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Abstract 

Based on a measurement of high momentum rf production in B decays, we determine the 
charmless inclusive B — > n'X nc branching fraction in the lab- frame momentum interval 2.0 < 
P rj i < 2.7 GeV/c. Using 9.7 x W®BB pairs collected at the T(4S I ) center-of-mass energy with the 
CLEO II and II. V detector configurations, we find B(B -» n'X nc ) = [4.6 ±1.1 ±0.4 ±0.5] x 10" 4 in 
the 2.0 < Ptf < 2.7 GeV/c momentum range, where the uncertainties are statistical, systematic, 
and from subtraction of background from B decays to charm, respectively. 
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The flavor changing neutral current decays b — > sg are forbidden at tree level in the 
Standard Model (SM), and hence at lowest order can only occur at the loop level. Interest 
in the charmless inclusive decay B — > rj'X (which we denote as B — > rj'X nc ) arises because it 
is expected to be dominated by the b — > sg transition followed by fragmentation of the gluon 
into if via QCD anomaly coupling arid formation of multiparticle states X by the s 

quark. The amplitude of these gluonic penguin decays may receive significant contributions 
from diagrams with virtual non-SM particles in the loop. 

CLEO previously reported an unexpectedly large rate for high momentum if production 
from B decays |8j. That result was based on 4.7 fb _1 of total luminosity, taken both on 
the T(4S) resonance and at a center-of-mass energy 60 MeV below the resonance, which is 
below the BB production threshold. These datasets are referred to as "on-resonance" and 
"off-resonance" respectively. CLEO found an inclusive B — > r]'X nc branching fraction of 
[6.2 ± 1.6(stat) ± 1.3(syst)l$;°(bkg)] x KT 4 for 2.0 < P rf < 2.7 GeV/c. 

In this paper we present a new measurement of if production from B decays in the mo- 
mentum range 2.0 < P n > < 2.7 GeV/c. Then, by subtracting background from B to charm 
decays, we use this result to determine B — > rfX nc . The result is based on 9.1 fbr 1 of 
on-resonance data and 4.4 fb _1 of off-resonance data. The analysis method is improved over 
CLEO's previous B — > r]'X nc analysis and now uses a combination of the "pseudorecon- 
struction" and neural-network/shape- variable approaches that have been used in the CLEO 
analyses of b — > (H, These strategies are used to isolate the signal and to suppress 
the contribution from continuum if production. We first search for B — > i]'X nc candidates 
that are consistent with one of the B meson multiparticle decays. We then estimate the 
background from B decays to charm via Monte Carlo technique. These new results in- 
clude the data used in the previous analysis and the results presented here supersede that 
measurement. 



The data used for this analysis were collected with the CLEO detector at the Cornell 
Electron Storage Ring, a symmetric e + e~ collider. The CLEO detector measures charged 
particles over 95% of 4ir steradians with a system of cylindrical drift chambers. (For 2/3 of 
the data used here, the innermost tracking chamber was a 3-layer silicon vertex detector.) 
Its Csl calorimeter covers 98% of Att. Charged particles are identified by specific ionization 
measurement (dE / dx) in the outermost drift chamber and by time-of-flight counters placed 
just beyond the tracking volume. Muons are identified by their ability to penetrate the iron 
yoke of the magnet. Electrons are identified by the ratio of their shower energy to track 
momentum, by track-cluster matching, and by their shower shape. 

We select events using a standard set of CLEO criteria for hadronic final states, and then 
search for candidate if mesons by reconstructing the if — > 7/7r + 7r~, 77 — > 77 mode, where 
the 77 forming the 77 candidate must have an invariant mass within 30 MeV of the nominal 
77 mass. We then form the mass difference between the reconstructed 77' and 77 masses to 
improve resolution. The mass difference M(it + tt~ 77) — M (77) must be within 50 MeV of 
the nominal mass difference AMp^g = M v > — M v = 410.5 MeV [2]. We restrict the 7/ 
momentum to P v > > 1.6 GeV/c. 

Using reconstruction we attempt to identify events in which a B decay produces a strange 
quark recoiling against an 7/. The reconstruction is done by forming combinations of a 
charged kaon or a K$ — ► 7r + 7r~, an if candidate, and n pions where n < 4 (at most one 
of these pions is allowed to be neutral); a total of eighteen decay modes and their charge 
conjugates are considered. For each B candidate we calculate the momentum P, energy E, 
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and beam-constrained mass M = yE 2 eam — P 2 . We then form the Xb °f the reconstruction: 

where cr E = 40 MeV and a M = 4 MeV. 

We only consider candidates with x B < 20, and if an event has more than one candidate 
with Xb < 20 we choose the candidate with the lowest Xb- This reconstructed X system 
may include strange, charm, or light quarks. Events with only light quarks in the recoil 
system (Xa( u )), however, will have slightly lower efficiency than the signal X s cases. Events 
with charm in the recoil system will have comparable efficiency to the X s cases, and we will 
later subtract these as background. In cases where all decay products of X are identified 
correctly (50% of the time) the X mass (M(X)) resolution varies from 20 to 30 MeV for 
low and high M(X) respectively. The M(X) resolution becomes 200 to 300 MeV if there 
are missing or extra particles in the X reconstruction. Given a reconstructed candidate, we 
use Xb an d I cos 9 a \ as variables for continuum suppression, where 6 tt is the angle between 
the thrust axis of the candidate B and the thrust axis of the rest of the event. 

We also use event shape variables and the presence or absence of a lepton (electron 
or muon) to further suppress the continuum background. Specifically, we use a neural 
network optimized on signal and continuum Monte Carlo to combine the following event 
shape variables into a single variable r s h ape : the normalized Fox- Wolfram second moment 
S± 1 , and the energies in 20° and 30° cones, parallel and antiparallel to the rj' 
direction. If the event contains a lepton then we also use the momentum of the lepton, Pi, 
and the angle between the lepton and the rj , 9i v '. 

We thus have two types of events: those with both a pseudoreconstruction and a lep- 
ton, for which we use r s h ape , Xb, I cos @tt\> Pi, an d 0i v '] and events with only a pseudore- 
construction, for which we use r s h ape , Xb, an d I cos 9 tt \- For each of these two cases the 
available variables are combined using a neural network that has been optimized using 
signal and continuum Monte Carlo samples. For each of the two resulting networks, we 
can maximize the statistical strength by converting the network output r into a weight, 
w(r) = s(r)/[s(r) + (1 + a)b{r)], where s(r) and b(r) are the expected yields for signal and 
for continuum background respectively, r is the net output, and a is the luminosity scale 
factor between on-resonance and off-resonance data samples (a ~ 2.0). 

The above choice of weights minimizes the expected statistical error on the B — > rj'X nc 
yield after off-resonance subtraction. To reduce the systematic error that arises from having 
some efficiency dependence on the reconstructed M(X) value, however, we adjust the weight 
event-by-event based on the measured value of M(X). The factor for this adjustment was 
obtained from signal Monte Carlo by fitting the measured efficiency dependence on M(X) to 
a straight line. By using these adjusted weights we slightly increase our expected statistical 
error, but decrease a systematic error. 

We obtain yields from both on- and off-resonance data by fitting the distributions of the 
mass difference M(7r + 7r~77) — M^j) — AMpdg- We use Gaussian and linear functions for 
the rj signal and combinatorial background respectively. In these fits, we constrain the mean 
of the Gaussian to be zero, and the width to be 4.0 MeV based on Monte Carlo simulation. 



S± is defined as the sum of the magnitudes of components of momenta perpendicular to the direction of 
the rj for all particles more than 45° from the rj axis, divided by the sum of the magnitudes of momenta 
of all particles other than the rj . 
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We scale the off-resonance yield by the on to off ratio of C/E^ m to account for luminosity 
and cross-section differences between the two datasets. We also scale particle momenta by 
the on/off energy ratio to account for the small energy difference between the on-resonance 
and off- resonance data 2 . We then find our B — > rj'X yield by subtracting this off- resonance 
yield from the on-resonance yield. The mass difference M(it + tt~ 77) — M(77) — AMp^c m 
the 7/ momentum range 2.0 < P v > < 2.7 GeV/c is shown in Figure ^ and yields from fits 
are tabulated in Table We estimate a systematic error of 3% from the uncertainty in the 
fitting procedure. 

To search for B — ► i]'X nc , we take as signal the momentum region 2.0 < P v i < 2.7 
GeV/c because it is above the kinematic limit for most B to charm decays. Some B to 
charm decays, however, will still enter this signal region. We estimate this background 
by using BB Monte Carlo events and by measuring the data yield in a control region of 
1.6 < P r] i < 1.9 GeV/c, chosen because here we expect B — > rj'X nc to be much smaller than 
B decays to charm. The estimated yields from B backgrounds are tabulated in Table 
The continuum-subtracted and combinatorial-background-subtracted M(X) distribution is 
shown in Figure 121 Cascade decays (e.g., B° — > Dt s ^X, — > rj'Y) and color allowed 
decays (b — > cW,W — > rj' + {nir)) have been simulated in the CLEO BB Monte Carlo. 
We do not use this Monte Carlo prediction directly; instead, we rescale the CLEO BB 
Monte Carlo yield in the signal region by the factor of 1.06 ± 0.22 to account for half of the 
difference between data and Monte Carlo in the control region. We include the statistical 
and systematic uncertainties from the control region in the error on the scale factor. We 
estimate the yield of color suppressed decays with a Monte Carlo sample that assumed 
the following branching fractions: B(B° -> r]'D° + rj'D *) = (3.1 ± 0.9) x 10~ 4 Q and 



B(B° -> r]'D **) = 1/2 B(B° -> r]'D° + rj'D *) = 1.55 x 10~ 4 . We take the error on 
B(B° -> ri'D **) to be 50% of itself. We use the modified ISGW model for the mass 
distribution of D **. 

We determine the efficiency — defined as weights per B — > v(X nc event generated in the 
momentum range 2.0 < P v / < 2.7 GeV/c — via Monte Carlo simulation. The hadronization 
of the quark pair into X s or Xm u ) is done by the JETSET Monte Carlo. We use a uniform X nc 
mass distribution for the multiparticle final states, and assume that the branching fraction 
for B — > rj'K is one eighth of the total branching fraction B — > i]'X s in this momentum 
range The detection efficiency is averaged over charged and neutral B mesons and 

corrects for unobserved modes with neutral kaons (K®, Kg — > 7r°7r°), for modes with a 
charged kaon and more than one ir° in the final state, and for final states with baryons; it 
also includes the product of branching fractions B(r)' — > i]n + 7i^) x B(r) — > 77) = 16.96%. 

To estimate the uncertainty in efficiency due to the choice of a uniform shape for the 
M(X nc ) distribution, we vary the shape of the M(X nc ) distribution for the multiparticle 
final states from uniform to linear with intercept at the K — ir threshold, keeping the fraction 
of B — > rj'K constant; this leads to a systematic error of 6.3%. The systematic error on the 
efficiency includes 3.3% uncertainty from the event modeling of B — > r]'X nc which includes 
uncertainties in event shape, hadronization, and other-S modeling, and 2.1% uncertainty 
from the detector performance which includes uncertainties in the detection efficiencies and 
resolutions for tracks and photons. 

To obtain the B — > i]'X nc branching fraction, we take the background-subtracted yield 
in the momentum range 2.0 < P v i < 2.7 GeV/c of 61.2 ± 13.9 ± 6.6 weights, where the first 



The on- and off-resonance datasets have center-of-mass energies of 10.58 and 10.52 GeV respectively. 
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error is statistical and the second is from uncertainty in the subtraction of the B to charm 
decay background. The efficiency is (6.81 ± 0.56) x 10~ 3 weights per event, where the error 
results from dependence of efficiency on the M(X nc ) distribution, from X nc hadronization, 
from other- B modeling, and from detector performance. Our sample contains 9.7 million 
BB pairs (±2%). We obtain B(B -> r]'X nc ) = [4.6 ± l.l(stat) ± 0.4(syst) ± 0.5(bkg)] x 10~ 4 
for 2.0 < Prf < 2.7 GeV/c. 

Also included in our measured branching fraction are components from b — > dg gluonic 
penguin decays and b — > u tree level decays. Within the SM, the contribution from these 
decays is expected to be of the order of a few per cent of the B — > r]'X nc rate. We determine 
that the ratio of efficiencies for the final states that consist only of d and u quarks to the 
final states that contain s quark is 0.79. Thus, the branching fraction we have measured is a 
weighted sum of branching fractions B{B — > rj'X nc ) = E[B — > r]'X s ) + 0.79 B(B — > r]'X d ^). 

In summary, we have updated the measurement of the charmless inclusive decay B — > 
r)'X nc for 2.0 < R n > < 2.7 GeV/c using the CLEO II and II. V data sets which contain 9.7 
million BB pairs. This result is in agreement with CLEO's previous inclusive measurement 
and hence confirms its unexpectedly large rate. CLEO has also previously measured [16] the 
exclusive decay B — > rj'K using the same dataset as the results we report here; that result 
is therefore not statistically independent of our new results. 

The large B — > 7/X n(j _rate can be understood if the b — > sg rate is larger than expected 
in the Standard Model [2j. A number of alternative explanations, however, have also been 
proposed. One of these is the QCD anomaly mechanism b — *> sg, g — > rj'g with the rj'gg form- 
factor being a slowly-falling function of the gluon Q 2 [lj. A slowly-falling r/'gg form- factor, 
however, is strongly disfavored by the recent CLEO measurement of the T(1S) — > rj'X 
spectrum [l7|; a rapidly-falling form- factor needed for consistency with the T(1S) — > r/'X 
result predicts a B — >• r/'X rate of 3 x 10~ 5 [18], which is an order of magnitude smaller 
than observed. Another possible explanation of the high rate for B — > g'X nc is that the rj' 
has a substantial intrinsic charm content 0. In this case the rf can be produced by the 
axial vector part of the b — > (cc)s process. This explanation is also disfavored by CLEO 
measurements. First, the B — > rj c K rate is not enhanced relative to the B — > J/ipK rate (l9| . 
Further, this explanation requires that the B — > r]'K* rate be roughly half the B —>■ rj'K 
rate 0,0, which is inconsistent with observations [l6l| . jiol] . 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. M. Selen thanks the Research Corporation, and A.H. 
Mahmood thanks the Texas Advanced Research Program. This work was supported by the 
National Science Foundation, and the U.S. Department of Energy. 
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rf Momentum Range (GeV/c) 


1.6 - 1.9 


2.0 - 2.7 


ON 

ax OFF 


251.8 ±25.9 
90.8 ± 15.3 


149.4 ± 11.4 
55.3 ±7.9 


ON-axOFF 


161.0 ±30.0 


94.1 ± 13.9 


Color suppressed decays 
B° -> r]'D° + r]'D°* 
B° -> r/ Z) ** 


0.0 ±0.1 
2.8 ±0.2 


17.4 ±5.0 
1.4 ±0.7 


CLEO Monte Carlo 
Scaled by control region 


140.6 ±11.0 


13.3 ±2.9 
14.1 ±4.3 


Sum of B backgrounds 


143.4 ± 11.0 


32.9 ±6.6 


ON— ax OFF— B backgrounds 


17.6 ± 30.0(stat) 


61.2 ± 13.9(stat) 



TABLE I: Yields(weights) from the fit in the signal (2.0 < P v > < 2.7 GeV/c) and control (1.6 < 
P v > < 1.9 GeV/c) regions. Given are yields on-resonance, scaled off-resonance, on minus scaled off, 
estimated background from color suppressed B decays, estimated background from cascade and 
color allowed B decays, sum of B backgrounds, and on minus scaled off minus B backgrounds. 
The error on the continuum-subtracted and S-backgrounds-subtracted yield is statistical only and 
does not include the error from the background subtraction. 
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FIG. 1: The distribution of M(ir + TT~jj) - M(jj) - AM PDG in the signal region 2.0 < P n > < 
2.7 GeV/c for (a) on-resonance data and (b) scaled off-resonance data. 
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FIG. 2: The continuum-subtracted and combinatorial-background-subtracted M(X) distribution 
(points with error bars) with estimated background from B decays to charm (histogram). The 
various contributions are cascade decays (light grey) and color suppressed decays B° — > rj'D°+r]'D * 
(dark grey) and B° — ► r]'D ** (hatched area). Ignoring the small smearing effects from the boosted 
B mesons (300 MeV/c), the mass ranges M(X) < 2.35 GeV and M(X) > 2.5 GeV correspond to 
momentum ranges in Table HI of > 2.0 GeV/c and P v > < 1.9GeV/c respectively. 
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